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ABSTRACT 


/ 


The instability of contrastreaming plasmas is investigated 
taking into account the Coulomb Collisions via the Fokker -Planck 
co-efficients in the Boltzmann equation. The dispersion relation 
is obtained on the assumption that the Coulomb collisions are 
weak and solved on the additional assumption that the phase 
velocity of the wave is much larger than the mean thermal velocity 
of the particles. It is found that while the temperature has the 
effect of increasing the maximum wavenumber (which for a cold 

plasma is equal to 0)‘k in units of , where 60^ is the electron 

plasma frequency and U is the streaming velocity) below which the 
plasma is unstable, the collisions have no effect on this wave- 
number. However, the growth rate of maximum instability decreases 
(compared to its value for a cold plasma) on taking into account 
the thermal motions but increases when the collisions are taken 
into account. 
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I. INTRODUCTION 


The instability of longitudinal electron oscillations 
in a plasma where we have two streams of electrons moving 
with equal and opposite velocities has been discussed exten- 
sively by several authors (l,2) when the thermal spread of the 
electrons as well as the collisions between the particles are 
ignored. It is well known that such a system is unstable for 
all wavenumbers x(= kv/ooi, , where & is the wavenumber of the 
perturbation, (J the speed of the stream and 6^the electron 

plasma frequency) which are less than a certain critical value 

\U 

=s (Sb) ' . Further the wavenumber for which the insta- 
bility is maximum is known to be = (3)^ • Jackson (3) has 

considered the problem of two stream instability when the 
velocities of the two streams are smeared out by an equal amount 

and finds that DC increases on taking into account the spread 

c 

of the distribution function of the electron streams. 

Tidman and Weiss ( 4 ) calculated the effect of collisions on 
the two stream instability assuming zero thermal spread of the 
electrons and the ions. The collisions were taken into account 
by a simple relaxation model in which the electrons are taken to 
be scattered into a local Maxwellian distribution in a time HZ . 
They found that collisions increased this instability. This 
result led them to suspect that the diffusion nature of the 
Coulomb thermalization of the electron streams is an essential 
feature of the problem and hence the more exact Fokker -Planck 
equation for Coulomb scattering should be used. It was later 
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pointed out by May (5 s ) that Tidman and Weiss did not carry out the 
calculations consistently as they neglected the equilibrium 
current density and the associated magnetic field which results 
due to the mass motion of the electrons. May considered two 
interpenetrating streams so that there is no current in equilibrium. 

He calculated the effect of collisions using the relaxation model 
of Gross and Krook (6) and found that the unstable mode ceases to 
grow after a time of the order of one collision time. However, he 
had taken into account only the electron -electron collisions 
claiming that the electron-ion collisions are negligible due to 
the large ion-electron mass ratio. We will show here that the 
electron-ion collisions are actually comparable with the electron- 
electron collisions for the problem at hand. 

Tidman ( 7 ) considered the effect of collisions onettie instability 
of a stream of particles moving with a velocity (with no thermal 
spread) through cold ions at rest. He took into account the effect 
of Coulomb collisions using the Fokker -Planck equation and found 
that as the ordered streaming energy of the electrons is thermalized 
due to diffusion in velocity space, there results an enchanced 
Landau damping which competes with the growth of instability. He 
did not take into account the equilibrium current density which 
results in his model. It will be shown from our more rigorous 
calculations that weak collisions cannot even cancel the effect of 
thermal motions much less quenching the instability. 

Comisar (8) and Buti and Jain (9) have studied the effect of weak 
Coulomb collisions on the damping of longitudinal and transverse 
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plasma oscillations respectively using the Fokker -Planck equation 
in the form given by Rosenbluth et al (lO). They find -that to the 

first order in the collision frequency l ) , the electron-ion 

C 

collisions dominate the electron-electron collisions in damping 
the plasma oscillations. 

Following Comisar and Buti and Jaw, we consider the effect of 
Coulomb collisions on the stability of contrastreaming hot plasmas 
in a systematic manner. We assume that the collisions are not too 
frequent so that we can make an expansion in powers of ) and 

we calculate the collisional effects to the first order in the 
collision frequency. In the absence of collisions, the dispersion 
relation, including the thermal effects, is exactly solved numerically 
to determine the region of instability and the growth rate of maximum 
instability. When collisions are taken into account along with the 
thermal effects, then we use an iterative procedure to solve the 
dispersion relation and to study the effect of collisions on the 
critical wavenumber and on the growth rate of maximum instability. 

It is, perhaps, in order to remark here that the Fokker -Planck 
equation is not strictly valid for an unstable plasma, particularly 
near the electron plasma frequency ( 11,12). However, the motivation 
for the present undertaking is to clarify the contributions of the 
frictional and the diffusion terms In the electron -e lectron and 
electron-ion collisions. We find that the contribution of the 
frictional term is dominant compared to that of the diffusion term 
in e 1 e c t r on - i on coll i r, i on ' ; v\ ,e r e a s i n e 1 o c t r on - e 1 o c l. r on collisions 

same order. Further, the contribution of the 


the two are of the 



electron-ion collisions is much larger than that of the electron- 


electron collisions. It may be reasonable to expect that the 
qualitative nature of these results will not be altered by using 
a more exact kinetic equation for an unstable plasma. 

II . GENERAL THEORY 

Let us consider two streams of hot unbounded plasmas in which 
the ions form only the neutralizing background but the electrons 
have non-relativistic streaming velocities. The equilibrium dis- 
tribution functions, normalized to unity, for the electrons and the 
ions are assumed to have the form: 

- (jl^ vf J ^ Av**. ^ 

and 


ii e 


a. ~ ^ 


(2) 


where Yx ~ /'**' and '=<3>./M ; ©. and denote the temperatures 

(in energy units) of the electrons and the ions respectively. The 
subscript a on the electron distribution function refers to either 
of the streams. The density of electrons on ions in each stream 
is taken to be N. The distribution functions given by (l) and 

(2) do not satisfy the Fokker -Planck equation: 


where F is the force per unit mass and 

w* r 

< 4 > = Nrz, (1 * %)£ ft y ' . 

oc=i,e C JLJ IV-YI 
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( 3 ) 


( 4 ) 


■* VIA w. 


r = ±i|i^y\ , 


where/l is the Debye length. In fact we find that 


i^?= ±2UM1 FaH(4fAv; a )+ 

a.t CTy fc cu»/vi) 3 L + J 


where 


i2(A=f u 4 e a ^u- , 


and the effective collision frequency is defined to be 


^ = C^/A) ^A 


In writing Eq. (7) we have taken the electron temperatures of each 

Si 3j 

stream to be the same and put V Q 

As we shall see later the effect of thermal motions enters the 

dispersion relation through the parameter {k.\)=A*' with A=%/v 
x - kU/ctij,. 

y^Now for all situations of physical interest ( R^Aq) must be much 
less than unity. Further we know that two-stream instability arises 
for values of X of order unity. Thus the condition implies 
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that/\ be less than unity for the unstable situations we are 
considering here. Since 


% 




( 10 ) 


we can write 




(n) 


or in order of magnitude 


j_ dLL = 3-3LfVo\ 
Oi^U dt 4 Oij, \(J J * 


(12) 


The quantity oLh ) is always much less than unity for 
physical situations of interest and we have seen that is 

also less than unity for the unstable situations under consideration 
here. Hence (J changes in time on a scale which is very much larger 
than a plasma period. On the other hand we know that in the absence 
of collisions the two stream instability has e-folding time of the 
order of a plasma period. Thus for the problem at hand, the dis- 
tribution function satisfies the Fokker -Planck equation to a good 
approximation and we may take {A in Eq. (l) to be constant. Moreover, 
the neglect oflhis term is justified on the ground that this is 
independent of X and thus would affect only the k~0 mode when we 
take the Fourier transform of the linearized equations in DC- space; 
and we are not primarily concerned with this mode here. In a similar 
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manner, it can be shown that the Fokker -Planck equation for the 
ions is also satisfied to an equally good approximation. 

For longitudinal oscillations the linearized Fokker -Planck 
equation for the motion of electrons has the form: 


ht 



A k E.V 




(13) 


where the right hand side represents the effect of collisions. 
Again the subscript a in Eq. (13) implies that this equation 
refers to either of the plasma streams, 1 or 2. The change in 
due to collisions results from collisions with the ions or 
electrons in either stream. Thus we can write 



where (X^S-rl or 2. The collision terms are given by 
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(17) 





&. = A ir£- f/ oi ( r) I y- y'I dv' 

'ox oVdV J 

(&f = awr A- f few ij£' , 

\/« ^vJ-nrrTT 


<A4£ = wr^ jh Cr;iy- r\ w . (ao) 

£ £ 

The expressions f or ^A^and ^4 can be written down from Eqs. (l9) 


id (20) 


^placing £>& by h . The electric field E in Eq. (13) 


is governed by the Poisson equation: 


v-e—^Z U(*.x.t)4!t. 


Equations (13) and (21 ) from the characteristic value problem 
which we shall solve by the Fourier -Laplace transform method. We 
thus define 

r°- C* r — ot-ik-x 

e (2: 


E (*). LtUoc E(z.t), (23; 

Um A j j vvu 

with (&lA^o £ or the integrals to be convergent. The Fourier -Laplace 
transforms of Eqs. (13) and (2l) are: 

*) - + i ' £•*£(.*. *) 

_ Hz. eja) ■ dfe =j^£k ) (24: 

on tv t'bt Jc*tl. 



and 




(25) 


where 9 , ( ) is the Fourier transform of the initial distribution 

function i .e . , 


CL P 

d k OC)=JJxe jjz.x.o). 


(26) 


In order to solve the characteristic value problem posed by Eqs. (24) 
and (25) > it is now convenient to take the Fourier -tran sform in 
velocity space; we thus define 




(27) 


(28) 


Equation (24) now takes the form 

In order to simplify the collision terms, we first observe that 
the Fourier-Laplace transform of Eq. (15) leads to 




_ 1 <r 

<2, 


On using the identies 

I I 




(29) 


and 


iK-ri 

V 


sit*. 


1 ■ci-Cv-x') , 
d f e ~ J. 


, , r,, , 


11 


we find that 




On using the foregoing results in Eq. (29) we obtain 



where 




Similarly 


Ysr-f^ ] 


gr-csr-? 

isr-qi 


(30b) 


(31) 


(32) 
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We shall now take & to be along the j? -axis and also assume 
that and further we set wrwr/L , where jLj represents the 

mean free path of the particles between collisions. Eq. (28) now 
reduces to 



Equation (36) can be readily integrated to give 

r ,, x /Af . ... / 

rij * 7 

- s' W V 


oO 


(37) 


r Q 


* z 1 ^ 

1 H 




oO 


+ 


where 


■/\ J 


sr'-(cc,a>, «;0. 

We shall now assume that the collisions are weak so that 
-kL»i. Under this assumption, to the lowest order, Eq. (37) 
gives 




&k!k 




(38) 
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-AGl'/h 

-e 


( 39 ) 


and 


(40) 


oo 

< % 

J? (4 ST.) = J drV'^L yof/k -grty-iofuA . 

06 * 

We now substitute the zero-order solution given by Eq, (38) into 

Eq. (37) to obtain: .r- eo / 

- / vl /- < \ i / v a # I / . v 


*{[kfat)+ ki&AjQJH) + Ki(^)QA^)} 

*[{K?(r!: q) + KJf, ?)] + K? (srl f) %(A, fjfl . 


From Eq. (25) we have 


E^) = ±Ti<LZ Fjk,o,*). te) 

On substituting for P in accordance with Eq. ( 41 ) into Eq. (42), 


we readily obtain 


where 


E (a) = 

* fa-j 


(43) 


+ 


J? 

rc. u : 

4 . n°° p (44) 

v)jQJ*4)+j3te , ?; q a (k. ?)fl 



( 45 ) 



We may note here that(^(k,s) depends only on the initial 
perturbation. If we consider only those perturbations for which 
£<k,s) is analytic in the complex s-plane, then for the Laplace 
inversion of Eq. ( 43 ), we have to consider only the zeros of^(k,s) 
which are given by 


(k,/>) =o . 


( 46 ) 


This gives us the desired dispersion relation. It is now convenient 
to define the effective collision frequencey as V 0 jj^ . The 

dispersion relation then reads 



III. EVALUATION OF THE INTEGRALS 


We first need to evaluate 




where ) . It is convenient to write this as 


where 


7t 


7i = 7 - % a ^ L = (V*. + 4 ^ • 

The integral which occurs in Eq. (49) can be expressed in terms 
of the error function and we find that 


(50) 


3, 


£ **K°£ri v >\- %)]?- - ^°i e +E i(7+!j> (51 } 

where = ot^ + / (*)*' <U 

E*f6*) = f <L%j , (52) 

DC 

which for large values of the argument has the asymptotic expansion 


E*l(«),dl X Jl- A+ -H- TT 

0 «w L J xT c2 3C 


£_ 

£ 


(53) 


JJX. L. «2 

On making the legitimate assumption that the phase velocity of 
the wave is much larger than the mean thermal speed i.e. x)/^o» 4, 
we obtain : 


3^Vq L + . . . 

Ccc-ku^r ' 


% i( - k ' o) ^ ( L kuj 3 -^ 


-h- 


(54) 


where we have put /6=-io)for the sake of convenience. 


We next wish to evaluate 



We may note that is independent of (J and as such of 1 §- . 

* o 

1,2. In order 


r ^ r ^ 

therefore, write 


We can , 

to compare the relative orders of the diffusion term and the frictional 
term, we split up Eq. (55) as follows: 


-*>(<%- yo/k. 




T 0 ' 

where i A y 


o* 


and 


represent the contribution of the 
frictional and the diffusion terms respectively to the electron-ion 
collisions. Further 




(57) 


and 


K u (0J2) r-lrf (58) 

w ^ jf* /o-i.-9i/+ ' 4 

On substituting for V oj^) and ) in accordan 

Eqs . (57) and (51) into Eq. ( 56 ) , we obtain for '• 


ce with 


r* _ e " 

If ait 2 vj 1 


oC 


Jds-e^y^ -± v^-±v?( oje,- 9 ;*j 




(59) 


It is now convenient to introduce the variable^ = <r'€ 2 —j5 


The integral 


■y 


can then be written as: 




(60) 
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Where now = 0 l A +v>(<r-h)M' k . Let us introduce a system 

of cylinderical polar co-ordinate-a Aa f* >i E <i- (60) then 


reduces to 


rSj 06 °o 


i =_ 


^ vr ^ y* 


JctcrJd^ **!>(- ~i v < 


a x {e /+ - s.d^CU)}S /f , 

where ^ is defined to be 

o“- f^c f a% 

S !f = i f 


2. a 

If we now put l'*h , this integral can be written as 

Kf < 

where '&■“( V a 4- vf )J 2. The integral takes the form 


r «2, o o OC oO 


I. = - e. 


0* fj) [1- Cfj] • «*> 


Again if we assume that f + »4 1, we can use the asymptotic 
expansion of error function and after some reductions we find 


A oo 


i = 
'f 


- A _ [iff *4 (v&iujHs ,< 

tvAJ^J £ ' 1 +M* Jar’ 


18 - 



where 


« 2 j. -•* 


1 - -f- 


3 Gig, ^°~ a - 


L M*- (M *-) 3 (M*f (p-f *) 7 J 

and we have put 

f ^ (£ ^ V ^° J )I V o ? 2* V °M> %r^^, (67) 

/*- - (yr\-h M)/nr, M . 

We now observe that is usually much larger than q, being of 
the order of the phase velocity to the mean thermal velocity. 

We can, therefore, carry out the expansions of the quantities 

p ^ 

appearing in the denominators in ^ . The resulting integrations 

are then elementary. We obtain after some straightforward cal- 
culations: 


+$ >}-¥-[*-% -ft j]> 

where we have retained terms up to ^ only, and 

A a = 1- yf/sL *%> , y- P - J- ( K°-+ ill). (69) 

On substituting the foregoing results in equation (65) and after 


some reductions we find 


lA' f 

J = -1 d£_ , 

'f -3 (w) /jl J 

/ 


-otkU\ 


S , 
3 / 


where 
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CO 

•?/ ' Sy J 1 ** 1 **Jj[(**® 


{i- 


A { +- y^oc. 1- (A, +- A a x.J a 

a ft f >, 3 At. 


and 


{<•- 


+ 


(A/ + 4a*)+ if 

IS- At 


I ~h Az 3c) a ^ (Ai-h Aa.X.y^' 


A '--k(*+ iU -)-£m- 


V^U n 


(72) 


The various integrals occuring in Eq, (71) can be evaluated in 
terms of the error function ( ^.) > where a )^ and 

the exponential integral Et/-X) . As is a quantity much 

larger than unity, being of order^one can use the asymptotic 
expansions of these functions. The calculations are lengthy 
but straightforward. The values of the various integrals are 


-L 


given in Appendix A. The result is (retaining terms up to order ) 

S a =J 3 . fi 


3 / 


2- AX 


~A.fi 


30-0. / j 

i Mi) 



• A, 1 

A? t 

Ai-J 

(73) 


~h 




3 At 




— C \ — £A? - 
A, L Af 



We now observe that J\ { can be written as 


4 -**(*-*£)' 


5> ; /' (74) 

where and jk /|) 1 1 . We can, therefore, 

again carry out the expansions of ^in powers of /k/>^ . 
Retaining terms up to the fifth power in 4/v we obtain : 
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T -LH- 3 ±IL>- !£ ) 

^3 > cx Vo a*J 


0» CL ft. ^ 

-IT r^. k a + 'iUa.f J-f* 1 - 3ik yL V) 3 ^ f>V 1 - 

a -*7 1 r ?ivj > vf 


^ - A £f/+*) - <76) 

o 

It is, perhaps, worthwhile to remark again that the only assumption 
which has been made to evaluate the integral i.^y* is that the phase 
velocity of the wave is much larger than the mean thermal velocity; 
and this is indeed a legitimate one for the problem at hand. 

We now assume that ^/M and <£<C l/^which is quite valid 
in all cases of physical interest; the second requirement, as a 
matter of fact, is a consequence of the first one. Under these 
assumptions, we have 


-£ a ~ -Lv* i> = r 
j * > r°- 


O- 1/ 
<x v o 


~J d y ■ 


The integrals have been evaluated numerically for suitable values 

of the parameter /\ ~ 7<r . We thus finally obtain: 


21 



Ol. 

1 . = 


A ib^_ [y (a , 3 kUo. , /g^V 7 

Or/* l/ 0 3 (fco- kU^) 3 i^ l K-kU^ (a-kUJ * j 

v- r ^ . /sfcVo' 1 


+ Z 




lik + li?Y° - fl+ nf/v.*)} (7?) 

**L>-*iy*'' • J5 

. • (I /» / / *W M / / ■*' I 


$ku 4- 


#■</; 


rWc^JaJ 1 - * 


(to-kUcfJ ? 


where we have put /0»-£fc>for the sake of convenience* Proceeding 


on similar lines, we obtain 


Ct 

7 . 

V - 


2£^' 


*£»« 


lA . / 3" 

'o 


(30) 


(jin )^K 3 (co- kU^L ^ ^ 

s l fa) — J 7 (co-hUt J) 

It is to be noted from Eqs. (79) and (80) that in ion-electron 
collisions it is the frictional term which is the dominant one; 
the diffusion term tends to be smaller by a factor or order 1 • 

Combining Eqs. (79) and (80) we obtain for the contribution of the 
electron -ion collisions: 

__ Aik 3 


h Saj. 3kyp- +_ h k vq z l 

(C0-^) x J 

_ V r ui , 7ku«__ , M&L (i + n£)} 


( 81 ) 


+-Z 


v: 


7 V 0 x (ix>-kuj 

Next we wish to evaluate 

o 




L 


If.J- 




(82) 


where 


K !f &}) - £ y—J, 




and 
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y 

The evaluation of JL^^yr*^ proceeds on similar lines as 
and we find: 

'«/ - sfes 


where 


^ = i(U»-U4)lv, , *--fa+iklL)/k K 


t-rt 






Combining Eqs. (85) and (86) , we obtain 


U -__! fp (A _ J° _ M \ , 7» (£ + id. 

’ 0n)H‘o?L 3L °t sio i ^ 


«• (89) 


In particular we find that for collisions between particles of 

the same stream, this reduces to 

ft/ ^ ^ v 

^ 3 4V' (9C 
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^ ( 91 ) 

/here 

l<3f (<*;%)= ^{-i v ?(r£*-yf-iii(<r-y t )}m 

ind 

^l/ r ‘ 2 ) = ~>f i “Va- yHtyr-y* )]■ (93) 

“ i . Hr 


into Eq. (91) , we 


I: 


7 £ 

On substituting for £»<% and K (°7£) ii 

f . - = _ * Jf**' e. * 1 J^[e^ t ~^oC^Cv+)] (9/ 


Ljl f Jf r e J t- l- 

fyl-i ^ K ^%- 

re now it is convenient to define 

y/ K ' u i)/ v ° 


X 

4 

where 

* t £ 

and 


■t'cr 


We now 

.a 


^ = a - + k % • 

first carry 


^ > .. ,. 

out the (^integration* The result is 


*<**£${%)} “/■{ %%%-± [r °7i- v ‘7l } : 


I ,_e 

3 */ 


(96) 


Again introducing a system of cylinderical co-ordinates in 


space and carrying out the angular integration we get 

( S, 


a, 0 c 




e /+ - <a u. 


£ ic»] 


(98) 


where 




<=o *• i 

-*V. 1 


(99) 


'Pi e T±ll 

f ° CvL-hyZ) 

After some elementary reductions, Eq. (99) can be put in the form 

where we have put . On substituting this expression 

for ^ into Eq. (98), we get 


( 100 ) 


cx- 

I.. 


/ r 


Vi 


- JC^cL 


'3if wtfjyjj 


F(* 


/ > 


( 101 ) 


/ 

where cL is defined in Eq. (87) and 


(102) 


R*> .[*_ (i/a ^ i (a £*{(?££)] 

X H^ ^~ xfh HC^) - *] . 

If we carry out the asymptotic expansions of the error functions, 


-r 

we obtain on retaining terms upto Cl : 

3*}- 


F«= 


— DC 


9“ 


a a x 


(* a +*X<*-*)* 


. (103) 
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Since Oc/a 1, we can further expand the denominators in 

this expression in powers of Ocj cc and obtain: 


FC*) 


rv 


OC 




at ore 




^-(-1 - b — +• 3 


3 


a \ 


a* 


4 - -3&L ( _l _ £-\ . jlL f- j_ 


a 




) 


) 


(104) 


^ w ~ ^ ~ -r* 

We now substitute this expression for p (3L) into 
given by Eq. (lOl). The resulting integrations are then elementary 


and we obtain : 

r* a. 


' 3 *f ’ 


fWi 


x 

t ( 
ex 


<\ a * 


3lS-- x -.3 - d [\ 

=2a| 2a/ 


+ lr 4 (0^%) 


(105) 


OJ » 1 H 


.) 


2 f A 

In a similar manner, we obtain: 


3U 




(#) ( 4 “ »;) ,r J • 


(106) 


On combining the foregoing expressions, we obtain: 

CL _ <x _ <V 


^ = 

_ 3 


1 


3-M 
r 


1 - 1 + ac 


AVfrS- _i-l' 1 (107) 


(j) (" l + c£ + <J; _£ i a *)} + ^7(a) {(‘k«A _a » 


o> at/ 


(<V ,( -a/ V 


3. 

X-C 

CX 


_ j 3 


aa^ a 


a. 

cv 


«3 


j 
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In particular for collisions in the same stream, this reduces to 


-r- 0* 

1 - 

3Cu 


30ir) /:2 ~v o € a} 

° OL 



(108) 


It is interesting to observe that for collisions in the same 
stream, the total contribution of the electron -e lectron collisions is 



£ _L 

• 3 * «■ t r^-vjaf ’ 

6 0 - 


(109) 


and compared with the electron-ion collisions, this is smaller by 
a factor Cl/ a . It turns out that the integral defined by 
Eq. (88) is about an order of magnitude smaller that the integral 
defined by Eq. (78). Hence the contribution of the electron- 
electron collisions between one stream and the other is an order 
of magnitude smaller than the electron-ion collisions. Further, 
in electron -e lectron collisions, the diffusion terms tend to be of 
the same order as the frictional terms. 

IV ANALYSIS OF THE DISPERSION RELATION 
On carrying out the summations indicated in Eq. (47) the dis- 
persion relation can be written as: 
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. We introduce 


where we shall now take 
the following dimensionless variables: 

a: * kU/o^ , ^ = o/ojp t A= K l M 


on using the foregoing relations, the dispersion relation 
be written as: 


1 . A Cft r£l + c. + C, * C, , , 


with 


C = CA* <? 

- siB . [f 

(jy/k ^t x ijZL J JL 1 0 

+ 1 

c T'^ 


c= 


and 


c, - - — ■■ # f- 4 — f V3xM/ ( 

+ 1, i- 4 * &, 


x<-j3 ’ ’ (yi-x.y*' 

__ a,f± T-Tl , 

l -2 ' ^ ^ ^ 


cr-*T 


( 111 ) 

can 

( 112 ) 

(113) 

(114) 


(115) 
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where 


Z 3 -7Z^%Z 7 )> (m) 

*"*/#• ^ /0 ^3 -V-ffe-tZ, -/Z)J, ( , 18 ) 

St ~ f ' ( 119 ) 



L - 3 T .±T _±T ^AT -h^T 
J* ‘ a 3 4 S £ ‘s' /) 7 ? ’ 


( 122 ) 



It is worth noticing here that /} ^ = (^/^) > a quantity which 

is always small compared with unity for all cases of physical 
interest. We shall now discuss some special cases of the dispersion 
relation (112). 


0 Of C a. 


c u 


In this case we obtain the well known cold case of two-stream 
instability which has been extensively discussed. The dispersion 
relation leads to the roots 




(125) 


The roots for 60 are thus either real or purely imaginary; the 

latter one gives rise to instability which occurs for all values of 

\L 

JC which are less than 3^“ &)' • It can be further seen that the 

%/ 

growth rate of instability is maximum for 3C. •= = (3) The 

p 

corresponding growth rate of maximum instability is I^J = o.a<t . 
Since the instability occurs only for the 'root with the negative sign 
in Eq. (78) , we shall henceforth consider only this root and denote 
it by the subscript 'O'. i.e. 


^ ~ li-oc x - (±-h^) 


% 


(126) 


We now wish to see how JC^and 3^. change when^^and Q'fco . 

b ) C x -Q, Cft = o, . » 

In this case the dispersion relation reduces to 


j _ + cfi .^ 


(127) 
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z 

This equation has been solved numerically for ^ for given 
values of 3L and ft . The results are given in Table 1 and are 
plotted in Fig. 1 to illustrate the effect of thermal motions on 
two stream instability. It is rather interesting to note that in 

this case the maximum wavenumber X ^ \ipto which the instability 

C 

occurs is larger thanXi.e. 


CO v 

= 4-4M-A. 


(128) 


It is to be further noted that while the region of instability is 
increased^ the growth rate of maximum instability is decreased, 

c) The effect of Collisions 

In order to compute the effect of collisions, we first observe 

that the terms Q and C are much smaller than unity so that we can 

Z /Z 

solve Eq. (ll2) by iteration. Thus we can substitute ^—2^ in 
where ^ is the solution of Eq. (127) and ^ is real for JC ^ 

Thus and (^become real in this region. Moreover ^ is positive 
for an unstable mode which we are considering here. Since and 
are proportional to^and at JC 9 we immediately conclude that 

the collisions cannot changeOJ and thus the region of instability. 

C 

However to study the influence of collisions on the growth rate of 
maximum instability, we have solved Eq. (112) for various values of 
A and B. The results are illustrated in Figures 2 and 3 • 

It is rather interesting to find that the growth rate of maximum 
instability, which is reduced by the thermal motions of the particles, 
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is enchanced by the collisions but the effect of thermal motions 
cannot be quenched by the Coulomb collisions. 

V. CONCLUSIONS 

In the absence of collisions, the region of instability of 
contrastreaming plasmas is increased by the thermal effects while 
the growth rate of maximum instability is decreased. However, on the 
assumption that the collisions are not too frequent l) 

one finds that they have only second order effects. The region of 
instability is not affected by the presence of collisions while the 
growth rate of maximum instability is increased. The binary collisions 
that we have considered here do not seem to quench the effects of the 
thermal motions of the plasma species. 
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APPENDIX 


We wish to evaluate integrals of the type 

1 *[<U ?±= Cil± ■ 

(A,****)* 


(A-l) 


These integrals can be readily evaluated in terras of the Error and the Exponential 
integrals defined by 


oo 

£*$(*) * Je^ cLy 


and 


Ei 


r" 

K-*)= - _e 

J SL 


y 3 


(A-2) 


(a-3) 


For large values of the argument, these functions have the asymptotic expansions (l 




and 




(A-4) 


(A-5) 


After some elementary reductions., we find: 

[l-A&etbfCi'j 


I = J. 

A A 1 

T i 


[l + £ e* Ei(rd* )1 » 


V i 


[1- *&e** E*{ CM] > 


<A 




K F f- , 

* U • v </# 




(A-6) 


(A-7) 


(A-8) 


(A-9) 
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and 


i. 


x 


* 0 ' A, A 


) Ax - 1 


where % = A</(2A X . These integrals have the asymptotic values: 



(A-10) 

(A-ll) 

(A-12) 

(A-13) 

(A-14) 





(A- 15 ) 
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TABLE 1 


A 

0.0 

0.001 

0.01 

0.05 

0.10 


B 


0.0 


0.000 

0.050 

0.075 

0.100 


0.000 

0.050 

0.075 

0.100 


0.000 

0.050 

0.075 

0.100 


0.000 

0.050 

0.075 

0.100 


x 

c 


1.4H 


1 . 416 
1.416 
1.416 
1.416 


1.435 

1.435 

1.435 

1.435 


1.516 

1.516 

1.516 

1.516 


1 .612 

l .612 
l .612 
l .612 


x„ 


0.866 


0.867 

0.867 

0.867 

0.867 


0.874 

0.874 

0.874 

0.874 


0.915 

0.916 

0.917 

0.918 


0.976 

0.982 

0.985 

0.988 


0.250 


0.249 

0.249 

0.249 

0.249 


0.245 

0.245 

0.245 

0.245 


0.227 

0.231 

0.232 

0.234 


0.214 

0.225 

0.230 

0.235 
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Figure Captions. 


Figure 1 . 


Figure 2. 


Figure 3. 


Title for 


Effect of thermal motions on x , x^ and y 2 for 

C " * 

A S 0.001, 0.01, 0.05 and 0.1. 

Variation of y 2 with x for B = 0 and 0.1 is illustrated 
for A - 0.05. 

Variation of y 2 with x for B = 0, 0.05 and 0.1 is 
illustrated for A a 0.1. 

Table 1 

The values of x , x, and y 2 for various values of 
c * * 

St 2 

A = tf/U 2 , and B = w /« . 

c p 
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+ 0.050 



- 0.300 











+ 0.050 



- 0.300 









+ 0.050 



- 0.300 








